MESH RATIOS FOR BEST-PACKING AND LIMITS OF MINIMAL 

ENERGY CONFIGURATIONS 



A. V. BONDARENKO*, D. P. HARDIN*, AND E. B. SAFF t 

Abstract. For Af-point best-packing configurations uj n on & compact metric 
space (A,p), we obtain estimates for the mesh-separation ratio j(ujn,A), which 
is the quotient of the covering radius of ujn relative to A and the minimum pair- 
wise distance between points in wjv. For best-packing configurations ujn that 
arise as limits of minimal Riesz s-energy configurations as s — > oo, we prove that 
-f(uJN, A) < 1 and this bound can be attained even for the sphere. In the particular 
case when N = 5 on S 2 with p the Euclidean metric, we prove our main result that 
among the infinitely many 5-point best-packing configurations there is a unique 
configuration, namely a square-base pyramid u)%, that is the limit (as s — > oo) of 
5-point s-energy minimizing configurations. Moreover, 7(^5, S 2 ) = 1. 

I. Introduction 

Let A be a compact infinite metric space with metric p : A x A — > [0, 00) and 
let un = {xi}f=\ C A denote a configuration of N > 2 points in A. We are chiefly 
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concerned with two 'quality' measures of ujn] namely, the separation distance of un 
defined by 

(1) S(u N ) = S p (u N ) := min p(x i ,x j ), 
and the mesh norm of un with respect to A defined by 

(2) r](u N ,A) =f] p (u N ,A) := max min p(y,Xi). 

y(zA l<i<N 

This quantity is also known as the fill radius or covering radius of u>n relative to 
A. The optimal values of these quantities are also of interest and we consider, for 
N > 2, the N -point best-packing distance on A given by 

(3) 5 N (A) = 5 P N (A) := max{5(u N ) : u N C A, #u N = iV}, 
and the N -point mesh norm of A given by 

(4) rj N (A) = V P N (A) := mm{r](u N , A) : u N C A, #w w = N}, 

where #S denotes the cardinality of set S. A configuration of iV points in A is 
called a best-packing configuration for A if 5(ojn) = 8n(A). 

In the theory of approximation and interpolation (for example, by splines or ra- 
dial basis functions (RBFs)), the separation distance is often associated with some 
measure of 'stability' of the approximation, while the mesh norm arises in the error 
of the approximation. In this context, the mesh- separation ratio (or mesh ratio) 

j(u N ,A) :=r)(u N ,A)/6(u N ), 

can be regarded as a 'condition number' for ujn relative to A. If {un}n=2 * s a 
sequence of iV-point configurations such that 7(wat, A) is uniformly bounded in N, 
then the sequence is said to be quasi-uniform on A. Quasi-uniform sequences of 
configurations are important for a number of methods involving RBF approximation 
and interpolation (see [9f [T3"l [T5| IT7]). 
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We remark that in some cases it is easy to obtain positive lower bounds for the 
mesh-separation ratio. For example, if A is connected, then 7(0;^, A) > 1/2. Fur- 
thermore, letting 

B(x,r) := {y £ A : p(y,x) < r} 

denote the closed ball in A with center x and radius r, then 7(0; at, A) > (3/2 for any 
iV-point configuration ojn C A whenever A and (3 £ (0, 1) have the property that 
for any r £ (0, diam(v4)] and any x £ A, the annulus -B(x, r) \ B(x, /3r) is nonempty. 
The diameter of A is defined by 

diam(A) := max{p(x, y) : x £ A, y £ A}. 

The outline of the paper is as follows. In Section 2 we present two simple but 
basic results concerning the mesh-separation ratio for best-packing configurations on 
general sets. In Section 3, we obtain lower bounds for this ratio for any best-packing 
configuration on the sphere in M. n and, in Section 4, we study the special case of 
minimal Riesz s-energy 5-point configurations on S 2 and determine their limiting 
best-packing configuration as s — > 00. Section 5 is devoted to a brief discussion of 
some special best-packing configurations on S n . 

2. Mesh-separation ratio for GENERAL SETS 
The following simple result is of the same spirit as that of Proposition 2.1 of [T2] . 

Theorem 1. Let (A, p) be a compact infinite metric space. Then, for each N > 2, 
there exists an N -point best-packing configuration ujn on A such that ^{ujn, A) < 1. 
In particular, this holds for any best-packing configuration ojn = {x\, . . . ,xn} having 
the minimal number of pairs of points {x{,Xj} such that p(xi,Xj) = 5n(A). 

Proof. Let be a best-packing configuration on A having the minimal number of 
unordered pairs of points {xi,Xj} such that p(xi,Xj) = 5jy(A). If r)(co^,A) > Sn(A), 
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then select a point a G A such that p(a,X{) > 5n(A) for i = 1, . . . ,N, and choose 
a point xt from some pair {x fc ,a^} such that p(xk,xi) = 6~n(A). Let u' N be the 
best-packing configuration obtained by replacing a in ojn by X£. Clearly, u' N has 
fewer unordered pairs of points {xi,Xj} such that p(xi,Xj) = 5n(A) than un- This 



On the other hand, there exist examples of compact metric spaces (A, p) for which 



as we now show. 

Example 1. Let A be the standard 1/3 Cantor set in [0,1] and let p be the Euclidean 
metric. For each iVeN, the set A is contained in the union of 2 N disjoint intervals 
of length 3~ N with endpoints = x^ < x 2 < . . . < x^ N+1 = 1 which belong to A. 
For any configuration of 2^ + 1 points in A, at least one of the intervals of length 3~ N 
must contain at least two points from the configuration showing that 5 2 n +1 (A) < 
3~ N . On the other hand, the configuration c^+i := {xf , ■ ■ ■ ,x^ N+1 = 2/3} is a 
best-packing configuration since S(u 2 n +1 ) = S 2 n + i(A) = 3~ N and has mesh norm 
t](u 2 n +1 ,A) = 1/3. Thus @ holds. 

Best-packing configurations arise as limits of minimum energy configurations as 
we now describe. For a configuration uj^ '■= {x\, . . . ,x^} C A of iV > 2 distinct 
points and s > 0, the Riesz s-energy of is defined by 



contradiction proves Theorem [TJ 



□ 



(5) 



limsup sup{7(co>7v ; A) \ 8(u N ) = 6~ N (A)} = oo, 



N^oo 




while the N -point Riesz s-energy of A is defined by 



(6) 



8 S (A, N) = €§(A, N) := M{E s (cj n ) : cu N C A, #lj n = N}. 
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An iV-point configuration un C A is said to s-energy minimizing if E s (un) = 
£s(A,N). 

Proposition 2 Q4|). Let (A, p) be an infinite compact metric space. For each fixed 
N > 2, 

U m £ s (A,N) 1/s = iATY 
Moreover, every cluster point as s — > oo of s-energy minimizing N -point configura- 
tions on A is an N -point best-packing configuration on A. 

The following theorem concerning the mesh-separation ratio of best-packing config- 
urations that arise as cluster points of s-energy minimizing configurations generalizes, 
simplifies, and improves Theorem 7 of [TT] . 

Theorem 3. For a fixed N > 2, let be a cluster point as s — >■ oo of a family of 
N -point s-energy minimizing configurations on a compact metric space (A,p). Then 
7(wn,A) < 1. 

The upper bound for j(u^,A) in this theorem can be attained even for the case 
when A is a sphere and p is the Euclidean metric. For N = 11 on S 2 , equality follows 
from the uniqueness result for best-packing of Boroczky [3J. For iV = 5 on S 2 , it 
follows from Theorem [7] in Section 4. 

Proof. Let iV > 2 be fixed and, for s > 0, let un, s be an iV-point s-energy minimizing 
configuration on A. Clearly, E s (un, s ) > S(un, s )^ s - This implies that there exists a 
point x s G con,s such that 

y£uN,s\{xs} 

lfr)(u) N)S ,A) > N 2/s 5(u N)S ), then E s (u)' Ns ) < E s (u) N>s ), where u' N>s := u N<s U{a}\{x s }, 
and a is a point of A such that pit, a) > rj(ojN tS ,A), for all t G un, s , which yields a 
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contradiction. Hence, 

(7) V("N, S ,A) < N 2 / s 8(u NjS ), 

and letting; s — > oo in and using Proposition [2j we obtain the statement of 
Theorem |3l □ 

3. Lower bounds for the mesh-separation ratio on the sphere 

In this section we derive some lower bounds for the mesh-separation ratio of a 
best-packing Appoint configuration on the unit sphere S n C IR n+1 with p the Eu- 
clidean metric. Let A n and G n be the sphere packing and covering constants in M n , 
respectively: 

(8) A n := lim N{5 N {U n )/2) n (3 n and 6 n := lim N VN (U n ) n Pn, 

where U n := [0, l] n denotes the unit cube in IR n and /3 n denotes the volume of the 
unit ball in W 1 (see, e.g. [7J [33]). First we prove the following asymptotic result for 
best-packing configurations on S n . 

Theorem 4. Let {un} denote a sequence of N -point best-packing configurations on 
S n . Then 

(9) ^/)>^(^ U + o(l) 1 N^oo. 



Proof. Since the collection of spherical caps with centers in the points of ojn of the 
radius t](ujn, S n ) covers S n , a standard projection argument implies 

(10) Np n (ri(u N ,S n )) n >e„Area(S") + o(l), N ^ oo. 
Similarly we have 

(11) Np n (^^J < A n Area(5") + o(l), N -> oo. 

Thus, we obtain ^ directly from (10) and (11). □ 
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It is interesting to investigate the asymptotic behavior of the constant on the right- 
hand side of ([9]) as n — > 00. The best known asymptotic upper bound for A n is the 
Kabatyanski-Levenshtein bound A n < 2 -0 - 599n+o ( n ) as n — > 00 and the best known 
lower bound for the covering constant is G ra > cn, where c is a positive absolute 
constant (cf. [3 pages 40 and 247]). Thus the inequality ^ implies the following: 
if n is large enough and N > C(n), then the inequality 

l{uJN, S n ) > (l/2)2 - 599 + o(l) > 0.757, n ->• oo, 

holds for an arbitrary best-packing configuration on S n . Further upper bounds 
for A n and lower bounds for B n can be found in [5] and [7j. In particular, it is known 
that for n = 2 the hexagonal lattice provides both A2 = vr /vl2, and G2 = 27r/v / 27. 
Hence 

l(ou N ,S 2 ) > -y= + o(l), iV^oo, 

for an arbitrary best packing configuration on S 2 . However, by special arguments 
working only for n = 2 we are able to improve this result to the following: 

Theorem 5. Let {lun} denote a sequence of N -point best-packing configurations on 
S 2 . Then 

(12) l(uj N ,S 2 )> 1 - + (l) = -J-= + (l), N^oo. 

2cos7T/5 1 + V 5 

Proof. It suffices to only consider sequences such that 7(co>tv, S 2 ) = 0(1) as iV — » 00. 
For a fixed N > 4, consider the Voronoi decomposition of S 2 generated by a; at, with 
Xi denoting the cell associated with xf, that is, 

Xi := {v E S 2 I \v — Xi\ = min \v — x\}. 

Euler's formula for convex polyhedra implies that there is a cell Xj having at most 5 
edges (each cell is a spherical polygon with edges consisting of arcs of great circles), 
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see [10]. Since 

B(x i} 5(u N )/2) n S 2 c X t c B(x i:V (uj N ,S 2 )), 

and rj(u)N, S 2 ) = 0(5(ujn)), it follows by a projection argument that there is at least 
one interior angle from Xj to consecutive vertices of Xj with angle 2tt/5 + o(l), and 
hence the distance from Xj to some vertex of Xj is at least 



2 cos 7r/5 



+ o(5(o;jv)), N -)■ oo. 



This yields (12). 



□ 



4. Limit of minimal energy for 5 points on S 2 




(a) (b) (c) 

FIGURE 1. 'Optimal' 5-point configurations on S 2 : (a) bipyramid BP, 
(b) optimal square-base pyramid SBP(l) , (c) optimal square-base 
pyramid SBP(16). 

It was observed in [15] from numerical experiments that 5-point minimum Riesz 
s-energy configurations on S 2 with the Euclidean metric appear to depend on s 
and to be of two general types: (i) the bipyramid (BP) consisting of 2 antipodal 
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points and 3 equally spaced points on the associated equator, and (ii) the square- 
base pyramid (SBP(s)) with one vertex at the north pole and 4 vertices of the same 
latitude depending on s and forming a square (see Figure [T]) . A comparison of the 
s-energy for the BP and the SBP(s) configurations is given in Figure [2] and suggests 
(as in [15] ) that BP is optimal for s < s* m 15.04808, while SBP(s) is optimal for 
s > s*. 

R. Schwartz [18J using a mathematically rigorous computer-aided solution proved 
(in a manuscript of 67 pages) that, for N = 5, BP is the unique minimizer of the 
Riesz s-energy for s — 1 and s — 2. (For the logarithmic energy, the optimality of BP 
is established in [8].) Currently there are no other values of s > for which a rigorous 
optimality proof is known. Regarding the stability of BP and SBP(s), in Figure [3] 
we plot the minimum eigenvalue of the Hessian of their s-energies. These graphs 
suggest that BP is not a local minimizing configuration for s > 21.148 (also observed 
by H. Cohn), while SBP(s) is not a local minimizing configuration for s < 13.5204. 

According to Proposition [2j every cluster point of s-energy minimizing configu- 
rations as s — > oo is a best-packing configuration. However, as is known, there are 
infinitely many non-isometric 5-point best-packing configurations on S 2 (see e.g. [2]). 

Proposition 6. ^(S" 2 ) = \/2 and all 5-point best-packing configurations on S 2 con- 
sist of two antipodal points (poles) and a triangle on the equator having all angles 
greater than or equal to tt/4. 

It appears from Figure [2] that the unique (up to isometry) cluster point of 5-point 
s-energy minimizing configurations is SBP(oo); that is, the square base pyramid with 
base on the equator. We shall next provide a rigorous proof that this is indeed the 
case. 
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Figure 2. The ratio of the s-energy of BP to the s-energy of SBP(s). 





(a) (b) 

Figure 3. The minimum eigenvalue of the Hessian for the s-energy 
of (a) the BP configuration and (b) the SBP(s) configuration. 

Theorem 7. Let Q' be a cluster point of a family of 5-point s-energy minimizing 
configurations on S 2 as s — > oo. Then Q' is isometric to 

(13) Q = SBP(oo) := {ei, -ei, e 2 , e 3 , -e 3 }, 

where ei = (1, 0, 0), e 2 = (0, 1, 0), and e 3 = (0, 0, 1). 
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It is perhaps surprising that this configuration has the maximum number of com- 
mon pairwise distances (eight) of length yj2 among all 5-point best-packings. 

We start the proof with an upper estimate for the minimum 5-point s-energy on 
S 2 . 

Lemma 8. 

Iimsup2 s / 2 £ S ( > S 2 ,5) < 8. 

s->oo 

Proof. For arbitrary < t < 1, we define the following 5-point configuration on S 2 : 

(14) Q t ■= {(±Vl-t 2 , -t, 0), (0, -t, ±Vl-t 2 ),e 2 }, 

which, for a suitable choice of t (depending on s), is a conjectured minimum energy 
configuration on S 2 for every s large enough. The s-energy of this configuration is 
given by 

E s {Q t ) := 4 • 2~ s (l - t 2 )- 3 ' 2 + 8 • 2- s ' 2 (l - t 2 )- 3 ' 2 + 8 • 2- s / 2 (l + t)~ s / 2 . 
Letting now t = s~ 2//3 , we obtain that 

lim(l -t 2 Y s/2 = 1 and lim (1 + t)~ s/2 = 0, 

and so 

Iimsup2 s/2 £ S (S 2 ,5) < lim 2 s/2 E s {Q t ) 
= lim (4 • 2~ s/2 (l - t 2 )- s/2 + 8(1 - t 2 )- s/2 + 8(1 + t)- s/2 ) = 8. 



□ 



We further need the following statement. 

Lemma 9. Let A, B, and M be fixed positive constants. Then 

f(x) := M{1 - Ax)' 8 + (1 + Bx)- S >M + min{l, AM/B} 
for every x G [0, l/A) and s > 0. 
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Proof. It is not difficult to see that / attains its minimum on [0, 1/A) at the point 
Xq = if B < AM and at the point 

(B/(M))^ S+1 ' - 1 
Xl ~ B + A(B/(AM)y/( s + l ) 

if B > AM. In the first case we have 

f(x)>f(0)=M + l, xe[0,l/A), s>0. 
In the second case, since 

*i < \ [(Bf(AM)) 1 ^ - 1] , 

we have 

f(x) > f(xi) > M + (1 + B Xi y s > M + (B/(AM))- s/(s+1) > M + AM/B 

for all x G [0, 1/A) and s > 0. Combining the results in both cases, we obtain the 
assertion of the lemma. □ 



Proof of Theorem [7l As we mentioned in Proposition [2] above, any cluster point 
of a family of s-energy minimizing configurations as s — > oo is a best-packing config- 
uration. Thus, by Proposition [6j it is sufficient to show that no 5-point configuration 
consisting of two opposite poles and an acute triangle on the equator (which we call 
an acute configuration) could be such a cluster point. We will prove this by contra- 
diction. For s large, consider a minimal s-energy configuration that is 'close' to a 
fixed acute configuration. We may assume that this minimal s-energy configuration 
Co>5(s) consists of three points 

A\ = Ai s = (ans, di2s, h), A 2 = A 2s = (a2is> &22s, h), A 3 = A 3s = (a^u, a^s, h), 
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where h = h s = o(l) as s — > oo, that are close to the vertices of a fixed acute triangle 
on the equator, and two points v4 4 = A± s and A 5 = A 5s that are close to (0, 0, 1) and 
(0, 0, — 1), respectively. Denote by 

3 3 

E x := E ls = J2\ A *- M~ S , and E 2 := E 2s = ^ \A b - A^ 3 . 

i=l i=X 

Clearly, the total s-energy E s (u 5 (s)) > 2Ei + 2E 2 . 

Let us first estimate E\ from below. Denote by O the point (0,0, h), by B the 
projection of A4 to the plane A%A 2 A^, and by x the length \0 — B\. Without lost of 
generality we may assume that B lies in the triangle OA 2 A^. Here we use the facts 
that x = x s = o(l) as s — > 00, and that AiA 2 A 3 is 'close' to a fixed acute triangle 
implying that O lies inside the triangle AiA 2 A 3 . Denote by a = a s , = (3 S , and 
7 = 7 S the angles A 2 OB, A 3 OB, A 2 OA\, respectively (see Figure [4]). 



Since 




i=l i=l 



we have, by the law of cosines and the fact that \B — A$\ = y/1 — x 2 — h, 
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E x = (2 - 2hVl -x 2 - 2xVl - h 2 



cos a 



-s/2 



+ (2 - 2h\fl - x 2 - 2xVl - h 2 cos f3)- s/2 

+ (2 - 2hVl - x 2 - 2xVl - h 2 cos(a + 7))~ s/2 . 

The crucial observation is the fact that a + (3 < t < tt, for some r that does 
not depend on s. Now monotonicity and convexity of the function t~ S//2 , t > 0, 
immediately imply 

(15) > 2 (2 - 2hVl — x 2 — x^l — /i 2 (cos a + cos /3) J 

+ (2 - 2Wl -x 2 + 2x)- s/2 

> 2 (2 - 2hVl - x 2 - xy/l ~ h 2 (l + cos t)) 

+ (2 - 2hVl - x 2 + 2x)- s/2 . 

From the facts that x = o(l), and h = o(l) as s — > 00 and the inequality 1 — x < 
— < 1, we get that 

E 1 >2(2-2h- #ixp /2 + (2 - 2h + 3x)~ s/2 , 

for some absolute constant Q\ > 0. Then, by Lemma [9j 

£1 > (2 + # 2 )(2-2/i)- s / 2 , 
for some absolute constant 9 2 > 0. Similarly we obtain 

E 2 > (2 + 2 )(2 + 2/i)- s / 2 , 

and so again applying the convexity of t~ s l 2 we finally deduce that, for s sufficiently 
large, 

(16) E s (u 5 (s)) > 2(E 1 + E 2 ) > (8 + A6 2 ) 2~ s / 2 . 
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On the other hand, from Lemma 8 we know that £ S (S 2 ,5) < (8 + o(l))2~ s / 2 . 



Therefore, by (16), an acute configuration cannot be a cluster point of minimal 



s-energy configurations as s — > oo. □ 

We can now obtain the dominant term in the asymptotic expansion for the minimal 
5-point s-energy. 

Theorem 10. We have 

lim 2 s ' 2 8 s (S 2 ,h) = 8. 

s— >oo 

Proof. By Lemma [8] it is enough to prove that 

(17) Iiminf2 s/2 £ S (S 2 ,5) > 8. 

For a fixed s > consider a minimal s-energy configuration Co>s(s) = {Ai, . . . , A^] : = 
{Ai s , . . . , A§ s }. By Theorem [7] we may assume that both distances \A2 — A^\ and 
|A 4 — A 5 \ have limit 2 as s — > oo. Observe that if the triangle AiA 2 A 3 is not acute, 
then 

|Ai - A 2 \- s + \A! - A 3 \~ s > \A\ - A 2 \- s + \A\ - A 3 \- s > 2 l - s ' 2 } 

where A[ is the midpoint of the circular arc (of length less than it) joining A 2 and 
A 3 and containing A\. A similar statement holds for triangle ^4x^44745. Therefore 
we may assume that at least one of the triangles Ax^^ or AiA^A 5 is acute since 
otherwise the desired lower bound for the s-energy follows. Without lost of generality, 
we assume that AiA 2 A 3 is acute. We adopt the same notation as in the proof of 
Theorem [7] and obtain a finer lower bound for E% and E 2 . 

There are three possible cases to consider, depending on the location of the pro- 
jection B of j4 4 onto the plane containing Ai, A 2 , and A 3 : (i) B is inside the sector 
A 2 OA 3 ; (ii) B is inside the sector A\OA 2 ] and (iii) B is inside the sector A\OA 3 . 
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Let us assume first that B is inside the sector A2OA3 as in Figure |4| From (15), 
we get 



E 1 > 2(2 - 2Wl - x 2 - xVl - h 2 (cos a + cos f3))~ s/2 > 2(2 - 2hy/l - x 2 ) 



-s/2 



In both other cases (ii) and (iii) we get the same inequality. Letting D denote the 
projection of A 5 onto the plane AxA 2 A 3 and setting y = \0 — D\, we similarly get 



E 2 >2(2 + 2h v /l-y 2 )- s / 2 . 

Thus, 

E s (u 5 (s)) > 4 [(2 - 2hVT^)- s/2 + (2 + 2h^/Y^Y s / 2 . 
Finally applying Lemma M to the last inequality and using the fact that x = o(l) 



and y = o(l) as s — )■ 00 we immediately obtain (17). □ 



5. Special best-packing configurations on S n 

In the case A = S n with n > 2 and Euclidean distance, there are best-packing 
configurations such that r](cuN, S n ) = 5(un) = \/2 for = n + 3, . . . , 2n + 1, 
yielding j(un, S n ) = 1 (see Theorem 6.2.1 [2]). For = 5 on S 2 , such a configuration 
is given by SBP(oo) defined in ( p~3j ). 

By the proof of Theorem [TJ we have t](un, A) > 5(lun) for some best-packing 
configuration un if and only if 5n(A) = Sn+i(A), which should be a very rare event, 
at least for A = S n . For S* 2 and N — 11 there exists a unique (up to isometry) 
best-packing configuration consisting of the regular icosahedron minus one of its 
vertices (see Hence, 

(18) Vn(S 2 ) = 5 n (S 2 ) and 7 (u> n , S 2 ) = 1. 

The unique best-packing configuration of 120 points on S* 3 is the 600-cell config- 
uration which has many other fascinating extremal properties, see [U [5]. Moreover, 
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in |19j . the numerical evidence is given that 

(19) 5 113 (S 3 ) = ... = 6 120 (S 3 ) = (V5 - l)/2. 



Assuming (19), we are able to construct a best-packing configuration of 113 points 
on the sphere with 77(0*113, S 3 ) > ^(wm). It consists of 600-cell without certain 7 
points which we describe below. 

In the 600-cell each point has 12 other points at the closest distance (a/5 — l)/2, 
and each pair of points at this distance has exactly 5 other points having the same 
distance to both points of the pair. So we will remove two points xi,x 2 , such that 
\x\ — x 2 \ = (a/5 — l)/2, and also 5 points yi, ■ ■ ■ 2/5, such that \xt — yj\ = (a/5 — l)/2, 
i = 1, 2, j = 1, . . . , 5. Recall that the second largest distance between points of the 
600-cell is 1. Thus, 



77(0113, S 3 ) > min 

xeojn-3 



Xi + x 2 

— X 



\X\ + x 2 \ 



3 1 ^ 1.2778 5(oii 3 ). 
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